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We investigate the thermodynamics and signatures of a polariton condensate over a range of 
densities, using a model of microcavity polaritons with internal structure. We determine a phase 
diagram for this system including fluctuation corrections to the mean-field theory. At low densities 
the condensation temperature, T c , behaves like that for point bosons. At higher densities, when T c 
approaches the Rabi splitting, T c deviates from the form for point bosons, and instead approaches 
the result of a BCS-like mean-field theory. This crossover occurs at densities much less than the 
Mott density. We show that current experiments are in a density range where the phase boundary 
is described by the BCS-like mean-field boundary. We investigate the influence of inhomogeneous 
broadening and detuning of excitons on the phase diagram. 

PACS numbers: 71.36.+c,42.55.Sa,03.75.Kk 



I. INTRODUCTION 

There have been many recent experiments with the 
aim of observing Bose condensation of polaritons in 
two-dimensional microcavities. Recent experimental 
progress includes nonlinear increase of the occupa- 
tion of the ground state-' 2 , sub-thermal second order 
coherence^, changes to the angular dispersion of polari- 
ton luminescenc o 3 i 4 i 5 increased population of low mo- 
mentum polaritons 3 -, and stimulated processes in reso- 
nantly pumped cavitiesSiiSiSiifl. Such experiments how- 
ever do not provide unambiguous evidence for the ob- 
servation of a polariton condensate as distinct from, e.g. 
a novel kind of laser. Theoretical predictions of the de- 
tailed properties of polariton condensation are therefore 
of considerable interest. In this paper, we consider the 
form of the phase boundary, and experimental signatures 
of condensation. We include the internal structure of po- 
laritons, and fcrmionic structure of excitons. 

The phase boundary at low densities may be described 
by a model of weakly interacting bosonsii* 1 ^* 1 ^. At higher 
densities, internal polariton structure becomes impor- 
tant, and the phase boundary approaches the result of a 
BCS-like mean-field theory. For a zero-dimensional cav- 
ity, the mean-field theory is correct at all densities, but in 
a two-dimensional cavity, fluctuations cause a crossover 
to the weakly interacting boson limit at low densities. 
This crossover is shown in the top panel of the phase 
diagram, figure The crossover scale is set by the wave- 
length of light, rather than average exciton separation, 
and so occurs at densities much less than the Mott den- 
sity. Figures 0] and also show the effect of detuning 
the exciton below the photon on the phase boundary. As 
in the mean-field case 1 ^, detuning can lead to a multi- 
valued phase boundary. However, with fluctuations this 
multi-valued structure occurs for smaller detunings than 
are required for the mean-field case. 

We also discuss a number of experimentally testable 
signatures of a coherent state of polaritons. These in- 
clude dramatic changes in the luminescence and absorp- 
tion spectra, as shown in figuresnand[21 and in the angu- 



lar distribution of radiation emitted from the condensate, 
as shown in figure [3J These signatures are a consequence 
of a coherent photon field, and of the presence of the 
Goldstone mode associated with symmetry breaking. In 
a previous publication 1 ^, some of these results were pre- 
sented in the absence of exciton-photon detuning. In this 
paper we present in full the calculation which led to those 
results, and extend our results to include detuning and 
inhomogeneous broadening of the excitons. 

We consider a model of localised excitons coupled 
to a continuum of radiation modes confined in a two- 
dimensional microcavity. This provides an extension of 
previous studies 1 ^ of the mean-field (zero-dimensional) 
system. A model of localised excitons is motivated by 
systems such as organic semiconductors 17-18 , quantum 
dot a 19 i 20 i 2 1 and disordered quantum wells 22 . In addition, 
the predictions of such a model are expected to be similar 
to those for a model of mobile excitons, since a typical 
exciton mass is several orders of magnitude larger than 
the typical photon mass. 

We study the behaviour of this model in thermal 
equilibrium. Although current experiments may remain 
far from equilibrium, there are several reasons to study 
the equilibrium behaviour. As the quality of micro- 
cavity fabrication, and particular the quality of mirrors 
improv e) 23 ! 24 , experiments can be expected to reach states 
closer to thermal equilibrium. Further, the nature of ex- 
perimental signatures for coherence close to equilibrium 
are expected to be similar to those in thermal equilib- 
rium. Finally, the equilibrium distribution may be con- 
sidered a limiting case of the non-equilibrium problem, 
when pumping and decay rates are taken to zero, so the 
equilibrium case is thus instructive in approaching the 
non-equilibrium problem. 

The Hamiltonian we study is similar to the Holland- 
Timmermans Hamiltonian 2 ^ 2 ^ studied in the context of 
Feshbach resonances in atomic gases. However, there are 
a number of important differences between the two mod- 
els; most notably the absence of a direct four-fermion 
interaction, and the bare fermion density of states. How 
these differences affect the mean-field results is discussed 
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in section UlI CI and their effect on the fluctuation spec- 
trum in sections 1111 Dl and 1111 El Because our system 
is two-dimensional, it is necessary to study fluctuations 
in the presence of a condensate. Such fluctuations have 
recently been studied by Ohashi and GriffinS! in the con- 
text of the Feshbach resonances. Despite differences in 
the models, we disagree with their use of partial deriva- 
tives of the free energy w.r.t. chemical potential, rather 
than full derivatives. As we show, their approach neglects 
terms of the same order as those that are included. 

The rest of this paper is organised as follows. In sec- 
tion [n] we introduce the Hamiltonian for our model, and 
discuss how it will be treated. Section IIIII reproduces 
the mean-field results presented elsewhere, and then dis- 
cusses the effective action, and resulting spectrum of fluc- 
tuations about it. From the fluctuation spectrum, a num- 
ber of experimental signatures are also discussed. 

In section llVl we explain how to calculate fluctuation 
corrections to the mean-field phase diagram. Some of 
the issues in this section are associated with the specific 
details of our model, or with fluctuations in two dimen- 
sions, but others are relevant to the calculation of fluc- 
tuation corrections in general. Section presents the 
results of including fluctuations, and discusses the effects 
of detuning and inhomogeneous broadening on the phase 
boundary. The discussion presented in section[V]does not 
require the details presented in section IIVI and readers 
not interested in the theoretical explication may move di- 
rectly to section IV1 Finally, in section IVII we summarise 
our conclusions. 

The calculation of the thermal Green's function for 
fluctuations in the condensed state includes terms pro- 
portional to Su- In the appendix, we discuss in general 
how such terms may arise in the thermal Green's func- 
tion, but vanish after analytic continuation and so are 
not present in the retarded Green's function. 



II. THE MODEL 

Our model describes localised two-level systems, cou- 
pled to a continuum of radiation modes in a two- 
dimensional microcavity. Considering only two levels 
describes a hard-core repulsion between excitons on a 
single site. We do not consider a static Coulomb in- 
teraction between different sites. The two-level systems 
may either be represented as fermions with an occupancy 
constraint, as will be described below, or as spins with 
magnitude |S| = 1/2. In the latter case, the generalised 
Dicke Hamiltonian^ is: 

j=nA 



H = E 2e ^ + E 

J = 1 fe=2iri/v / A 



+ ^=E (f^^Sj + e-^^lSj) ■ (1) 

j,k 



Here A — > oo is the quantisation area and n the areal den- 
sity of two-level systems, i.e. sites where an exciton may 



exist. Without inhomogeneous broadening, the energy of 
a bound exciton is 2e = %uj$ — A, defining the detuning 
A between the exciton and the photon. When later an 
inhomogeneously broadened band of exciton energies is 
introduced, A will represent the detuning between pho- 
ton and centre of the exciton band. The photon disper- 
sion, for photons in an ideal 2D cavity of width w, and 
relative permittivity e r is 



hw k = %-^=\ k 2 + ( — 
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2m ' 



(2) 



so the photon mass is m = (Ji^fe^ / 'c)(2ir j 'w). 

The coupling constant, g, written in the dipole gauge 

is, 
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2e e r 



(3) 



where d a b is the dipole matrix element. For small photon 
wavevectors (w.r.t. 1/w), it is justified to neglect the k 
dependence of g, i.e. to take u)k = ^o- The factor of 
1/y/w is due to the quantisation volume for the electric 
field. 

The grand canonical ensemble, H = H — /j,N, allows 
the calculation of equilibrium for a fixed total number of 
excitations, N, given by; 



N 



n A 

E 



57 



E 



^k- 



(4) 



We therefore define U&u = hu>k — H and e = e — fi/2. 
Expressing all energies in terms of the scale of the Rabi 
splitting, g\/n, and all lengths via the two-level system 
density n, there remain only two dimensionless parame- 
ters that control the system, the detuning A* = A/g^/n 
and the photon mass m* = mg/h 2 ^/n. Typical val- 
ues for current experiments^ are g\/n « lOmeV, m « 
I0 _5 TOeioctron, and taking n w l/a| ohr ~ 10 12 cm~ 2 leads 
to an estimate of m* w 10 -3 . 

This model is similar to that studied by Hepp and 
LiebSi. They considered the case n = 0, i.e. without a 
bath to fix the total number of excitations. The tran- 
sition studied by Hepp and Lieb was later shown by 
Rzazewski et al. 30 to be an artefact of neglecting A 2 
terms in the coupling to matter. The transition in our 
model, when \i ^ does not suffer the same fate^: The 
sum rule of Rzazewski et al. which prevents condensation 
no longer holds when fi =/= 0. 

In order to integrate over the two-level systems, it is 
convenient to represent them as fermions, S z — \(tfib — 
a^a) and S + — b^a. For each site there then exist four 
states, the two singly occupied states, |0), b^ |0), and 
the unphysical states |0) and a^b' |0). 

Following Popov and FedotoA^ the sum over states 
may be restricted to the physical states by inserting a 
phase factor e*^/ 2 ^ 6 b+ata ). Since the Hamiltonian has 
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identical expectations for the two unphysical states, this 
phase factor causes the contribution of zero occupied and 
doubly occupied sites to cancel, so the partition sum in- 
cludes only physical states. Such a phase factor may then 
be incorporated as a shift of the Matsubara frequencies 
for the fermion fields, using instead u) n = (n+ 3/4)2ttT. 
Thus, from here we shall describe the two-level systems 
as fermions. 



A. High energy properties, ultraviolet divergences 

The Hamiltonian Q is a low energy effective theory, 
and will fail at high energies. This theory is not renor- 
malisable; there exist infinitely many divergent one par- 
ticle irreducible diagrams, and so would need infinitely 
many renormalisation conditions. These divergent dia- 
grams lead to divergent expressions for the free energy 
and density. 

To treat this correctly, it would be necessary to re- 
store high energy degrees of freedom, which lead to a 
renormalisable theory. Integrating over such high energy 
degrees of freedom will recover, for low energies, the the- 
ory of eq. (JJJ. One may then calculate the free energy 
for the full theory, with appropriate counter terms. The 
low energy contributions will be the same as before, but 
the high energy parts differ, however such high energy 
parts are not relevant at low temperatures. As we are 
interested only in the low energy properties of this the- 
ory, we will introduce a cut off K m . The coupling g is 
assumed to be zero between excitons and those photons 
with k > K m . 

A number of candidates for this cutoff exist; the reflec- 
tivity bandwidth of the cavity mirrors, the Bohr radius 
of an exciton (where the dipole approximation fails), and 
the momentum at which photon energy is comparable 
to higher energy excitonic states (for which the two-level 
approximation fails). Which of these scales becomes rel- 
evant first depends on the exact system, however changes 
in K m will lead only to logarithmic errors in the density. 



III. MEAN FIELD AND FLUCTUATION 
SPECTRUM 

A. Summary of Mean Field results 



ipo, satisfies the equation: 



We first briefly present the results of ref s . IT 4l 1 6l for the 
mean-field analysis. Integrating over the fermion fields 
yields an effective action for photons: 

JO u 



M~ 



d T +l 



3s St 1 
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d T — e 



We then proceed by minimising S[ip] for a static uniform 
if) and expanding around this minimum. The minimum, 



hui tp = g 2 n 



tanh(/3£) 
2E 



■00, 



■r 



. 1-00 1 



.4 



(6) 



which describe a mean-field condensate of coupled co- 
herent photons and exciton polarisation. The mean-field 
expectation of the density is given by 



Pm.f. 



A 



1 



E 



taah(PE) 



(7) 



Note that the photon field acquires an extensive occupa- 
tion, we may define the intensive quantity po = l^o) 2 /^ 
the density of photons in the condensate. This corre- 
sponds to an electric field strength of yjTiLOQp$/2EQ. 

For an inhomogeneously broadened band of exciton en- 
ergies, eq. (0 should be averaged over exciton energies. 
In this case condensation will introduce a gap in the ex- 
citation spectrum of single fermions^: Whereas when 
uncondensed there may be single particle excitations of 
energy arbitrarily close to the chemical potential, now 
the smallest excitation energy is 2g\^)\j \J~A. The exis- 
tence of this gap is reflected by features of the collective 
mode spectrum discussed below. 



B. Connection to Dicke Superradiance 

It is interesting to compare this mean-field condensed 
state, described by eq. © , to the superradiant state orig- 
inally considered by Dicke2£. Dicke considered a Hamil- 
tonian similar to eq. but with a single photon mode. 
Constructing eigenstates \L,m) of the modulus and z 
component of the total spin, S = V. Sj, Dicke showed 
that for N particles, the state |JV/2,0) has the highest 
radiation rate. 

Taking the state described by the mean-field conden- 
sate, eq. ©J and considering the limit as ipo ~~ > 00 an d 
T — > 0, the equilibrium state may be written as 



IV'cond. 



4ll(lt> 




(8) 
(9) 



p=0 



i.e. a binomial distribution of angular momentum states. 

As N tends to infinity, this becomes a Gaussian cen- 
tred on the Dicke super-radiant state, with a width that 
scales like y/~N. However, it should be noted that the 
above state represents only the exciton part; this should 
be multiplied by a photon coherent state. In the self 
consistent state, there exists a sum of terms with differ- 
ent divisions of excitation number between the photons 
and excitons. These different divisions have a fixed phase 
relationship. Such a statement would remain true even 
if projected onto a state of fixed total excitation num- 
ber. This is different from the Dicke superradiant state, 
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which has no photon part. In the Dicke state, the only 
important coherence is that between the different ways 
of distributing excitations between the two- level systems. 



C. Comparison of Mean-Field results with other 
boson-fermion systems 

The mean- field equations ©, Q have a form that 
is common to a wide range of fermion systems. How- 
ever, the form of the density of states, and the nature of 
fermion interactions can significantly alter the form of the 
mean-field phase boundary. It is of interest to compare 
our system to other systems in which BCS-BEC crossover 
has been considered, and to note that even at the mean- 
field level, important differences emerge. To this end, we 
compare the mean-field equations for our modified Dicke 
model to the Holland-Timmermans Hamiltonian 25 - 26 27 , 
and to BCS superconductivity^. 

For a gas of fermionic atoms, or BCS, the density 
of states is non-zero for all energies greater than zero, 
whereas for our localised excitons, without inhomogc- 
neous broadening, it is a S function. One immediate 
consequence is a difference of interpretation of the mean- 
field equations; the number equation and the self consis- 
tent condition for the anomalous Green's function. For 
a weakly interacting fermion system the number equa- 
tion alone fixes the chemical potential^, and the self- 
consistent condition can then be solved to find the criti- 
cal temperature. For our localised fermions, the chemical 
potential lies below the band of fermions, and so the den- 
sity is controlled by the tail of the Fermi distribution, so 
temperature and chemical potential cannot be so neatly 
separated. This can remain true even in the presence of 
inhomogeneous broadening; the majority of density may 
come from regions of large density of states in the tail of 
the Fermi distribution. 

Such differences are also reflected in the density de- 
pendence of the mean-field transition. The absence of 



J 



a direct four-fermion interaction means that the effec- 
tive interaction strength is entirely due to photon medi- 
ated interactions. Since our model has the photons at 
chemical equilibrium with the excitons, this effective in- 
teraction strength depends strongly upon the chemical 
potential. 

For BCS superconductors, and the BCS limit for 
weakly interacting fermionic atoms, the dependence of 
critical temperature upon density is due to the chang- 
ing density of states, which appears in the self consis- 
tent equation as a pre-factor of the logarithm in the BCS 
equation l/g = p s hi(u) d/T c ). In contrast, the density de- 
pendence of T c for the Dicke model is due to the chang- 
ing coupling strength and occupation of two-level sys- 
tems with changing chemical potential. Such a change 
of coupling strength with chemical potential also occurs 
in the Holland-Timmermans model, near Feshbach reso- 
nance, where the boson mediated interaction is not dom- 
inated by the direct four-fermion term. Even with inho- 
mogeneous broadening of energies, unless the chemical 
potential remains fixed at low densities, there will be a 
strong density dependence of T c as the effective inter- 
action strength changes. The resultant phase boundary 
with inhomogeneous broadening, at low densities, is given 
in section IV"Cl 



D. Effective action for fluctuations 

Including fluctuations about the saddle point, tp = 
tpa + Sip, one may write the two-level system inverse 
Green's function as M -1 = M^ 1 + SM^ 1 , where 
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thus eq. (JSJ can be expanded to quadratic order as 



S = 
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(11) 
(12) 



The matrix between the photon fields can be identified 
as the inverse of the fluctuation photon Green's function, 
G 1 , where the exciton contribution to the quadratic 
term (with v=(n + 3/4)2ttT) is 



(iv + e) (if + iuj — e) 
[v 2 + E 2 ) {(v + uj) 2 + E 2 ) 
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where the sum over sites assumed no inhomogeneous 
broadening, and w is a bosonic Matsubara frequency, 
2imT '. The mean-field action, S[ipo], is given by: 



S[ip ] = Uoj k \ip \ 



N 
1 



In [cosh (j3 E)] . (16) 



The terms aS^, occur when the sum over fermionic fre- 
quencies in eqs. (TT3|l l(Tf|> have second order poles. These 
terms must be included in the thermodynamic Green's 
function at u> = 0. However, they do not survive ana- 
lytic continuation, and so do not appear in the retarded 
Green's function or in the excitation spectrum. This is 
discussed in Appendix El 

Even considering inhomogeneous broadening of exci- 
ton energies, such terms remain as 5 U , rather than some 
broadened peak. This can be understood by considering 
which transitions contribute to the excitons' response to 
a photon, i.e. between which exciton states there is a 
matrix element due to the photon. If uncondensed, the 
photon couples to transitions between the exciton's two 
energy states, ±e. In the presence of a coherent field, 
these energy states mix. The photon then couples both 
to transitions E — » —E and also to the degenerate tran- 
sition E — > E. Since this transition is between the two 
levels on a single site, inhomogeneous broadening does 
not soften the 5 U term. 

This conclusion differs for models with transitions be- 
tween two bands of fermion states. If transitions are al- 
lowed between any pair of lower and upper band states, 
the degenerate transition above is replaced by a contin- 
uum of intra-band transitions. In our model all intra- 
band excitations are of zero energy. If there exists a 
range of low energy intra-band transitions, these allow 
the Goldstone mode to decay, giving rise to Landau 
damping2L2i. For our model, as there is no continuum 
of modes, no such damping occurs. 

In order to consider fluctuations for an inhomoge- 
neously broadened system of excitons, a correct treat- 
ment requires calculating for a given realisation of dis- 
order, and then averaging the final results. Because the 
position of an exciton matters in its coupling to light, 
it would be necessary to include averaging over disorder 
in both energy and position of excitons. However, for 
low energy modes, we may make an approximation, and 
average the expressions for K\, K 2 over exciton energies. 

This approximation is equivalent to the assumption 
that the energies and positions of excitons sampled by 
photons of different momenta are independent and un- 
corrected. Such an approximation evidently cannot hold 
for high momenta, as otherwise the number of random 
variables would become greater than the number of ex- 
citons. This approximation will also necessarily neglect 
scattering between polariton momenta states. However, 
such effects involve high energy states (since they require 
momenta on the order of the inverse exciton spacing), and 
can be neglected in discussing the low energy behaviour. 



E. Fluctuation Spectrum 

Inverting the matrix in the effective action for fluctua- 
tions, eq. (|12f> . one finds the fluctuation Green's function. 
The location of the poles of this Green's function give 
the spectrum of those excitations which can be created 
by injecting photons, measured relative to the chemical 
potential. 

These poles come from the denominator: 

det^ 1 ) = \iu + hu k + Ki(uj)\ 2 -\K 2 (u)\ 2 



(lo 2 + AE 2 ) 



(17) 



where, as discussed above, we have ignored the 6 U terms, 
and assumed no inhomogeneous broadening of excitons 
In the condensed state the poles are: 



f 2 
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l -{A(k)±^A(k) 2 -B(k)) 1 (18) 



where 



A(k) = 4E 2 + (hu k ) 2 +4ihQ Q , (19) 
fi 2 k 2 

B(k) = 16 (E 2 m k ~ e 2 TuI> ) . (20) 

2m 

In the normal state this simplifies further, K 2 is zero, 
and eq. I|17(l is replaced by 



\iuj + huj k + Ki(uo)Y 



(iu> + E + )(iuj + £L) 



(iu> + 2e) 



(21) 



There are two poles, the upper and lower polariton; 



E± = - ((h£u k + 2e) ± yj(huj k - 2e) 2 + Ag 2 n tanh(/?e) 

(22) 

The polariton dispersion (1221) from localised excitons has 
the same structure as from propagating excitons, since 
the photon dispersion dominates. The spectra, both con- 
densed and uncondensed, are shown in figure ^ 

The difference between condensed and uncondensed 
spectra is dramatic: two new poles appear. These arise 
because the off diagonal terms in eq. I|12f> mix photon cre- 
ation and annihilation operators. Such a spectrum may 
be observed in polariton condensation experimentally, as 
one may probe the response to inserting a real photon, 
and observing its emission at a later time. In the presence 
of a condensate, the polariton is not a quasi-particlc. Cre- 
ating a photon corresponds to a superposition of creating 
and destroying quasi-particles. At non-zero temperature, 
a population of quasi-particles exists, so processes where 
a quasi-particle is destroyed is possible. This means that 
a process in which a photon, with energy less than the 
chemical potential, is added to the system is possible (i.e. 
^absorb as defined in eq. (|A7|) ~). However, if the system is 
experimentally probed with photons of energy less than 
the chemical potential, what will be observed is gain, 
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FIG. 1: Excitation spectra in the condensed (grey) and uncon- 
densed states, superimposed on the mean-field phase diagram, 
to show choice of density and temperature. The upper figure 
is for the photon and exciton resonant, and the two spectra 
are for T = 2<?y / n and fi/g^/n = —1.4 and —0.24. The lower 
has the exciton detuned by 3g\/n below the photon, and the 
spectra are for T — 0.1g^/n and fi/g^/n = —3.29, —3.01, 
— 2.54 and —0.37. The photon mass is m* = 0.01. Tempera- 
tures and energies plotted in units of g\/n, densities in units 
of n, and wavevectors in units of y/n. 



since the absorption of photons is more than cancelled 
by spontaneous and stimulated emission. 

At small momentum, £1 corresponds to phase fluctu- 
ations of the condensate, i.e. it is the Goldstone mode, 
and has the form £1 = ±hck, with, 

/ 1 { 4hu, g*n \ f\^\ _ [T~^ 
C -^2rn{ 6(0) 2 ){ N ) ~ V 2m n ' 1 ' 

The second expression is, for comparison, the form of 
the Bogoliubov mode in a dilute Bose gas, of interac- 
tion strength A. As tpo increases, the phase velocity first 
increases, then decreases. The decrease is due to the sat- 
uration of the effective exciton-photon interaction. 

The leading order corrections £1 = ±ck + uk 2 are 
of interest for considering Beliaev decay of phonons^. 
If a < 0, kinematic constraints prevent the decay of 
phonons. For the Bogoliubov spectrum in a dilute Bose 
gas, a = 0, and the cubic term becomes relevant. Here, 
both signs are possible; according to whether the spec- 
trum crosses over to the quadratic lower polariton disper- 
sion before this crosses over to a fiat exciton dispersion. 
In most cases, a > and £i will have a point of inflec- 
tion, but if the exciton is detuned below the photon, then 
for certain densities, a < 0, and the curvature is always 
negative. 

The modes may also be compared to the Cooperon 
modes in BCS theory. At small momenta, although the 
mode £i becomes a pure phase fluctuation, the other 
mode £2 is not an amplitude fluctuation. To see why this 
is so, it is helpful to rewrite the action in equation H12|) 
in terms of the Fourier components of the transverse and 
longitudinal fluctuations of the photon field. These com- 
ponents, at quadratic order, are equivalent to the phase 
and amplitude components, and are given by: 

Slpu.k = i>L{u,k) +iip T (to,k), 

Sipoj.k = ipL(-v,-k)-iipr(-v,-k). (24) 



I 

In terms of these new variables, the action may be written: 

iuj,k 

1 



Due to the off-diagonal components, the eigenstates 
are mixed amplitude and phase modes. This mixing van- 
ishes only where u) is small, which means that the low- 
est energy parts of the Goldstone mode are purely phase 
fluctuations. Since the amplitude mode at k = has a 
non-zero frequency, it will mix with the phase mode. 

The off diagonal terms come from two sources. The 
dynamic photon field leads to the term u>. The fermion 



mediated term ^s{K\) will be non zero if the density of 
states is asymmetric about the chemical potential. For 
the Dicke model, both terms contribute to mixing since: 

2 tagh(/lg) ( euj \ 

Unless e = 0, the density of states is not symmetric about 
the chemical potential, and so this term is non-zero. 
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For BCS superconductivity the pairing field is not dy- 
namic, so there is no off diagonal lo term. However there 
can still be mixing due to Sy(i^i), which is given by: 



iv.q 



C q+Q 



in which v is a fermionic Matsubara frequency, (2n + 
1)ttT, and A the superconducting gap. Note that this 
coupling now depends on the momentum transfer Q as 
well as energy lo. If the density of states is symmetric, 
e.g. e q = «f<Z, then at Q = 0, ^s{K\) will be zero, and 
as the boson field has no dynamics, the modes will then 
entirely decouple. In real superconductors, symmetry of 
the density of states about the chemical potential is only 
approximate, and so some degree of mixing will occur. 



F. Luminescence spectrum 

Although four poles exist, they may have very different 
spectral weights. At high momentum, the weights of all 
poles except the highest vanish, and the remaining pole 
follows the bare photon dispersion. 

Such effects can be seen more clearly by plotting the 
incoherent luminescence spectrum. As discussed in ap- 
pendix ^ this can be found from the Green's function 
as 



Pe 



emit fa) = 2n B (x)^ [Goo(i^ = x + irj)] . (28) 

It is easier to observe how the spectral weight associated 
with poles changes after adding inhomogeneous broad- 
ening. Figure |3 plots the luminescence spectrum for the 
same parameters as in figure ^ but with inhomogeneous 
broadening. This will broaden the poles, except for the 
Goldstone mode (labelled (d) in figure , as discussed 
above. The distribution of energies used is, for numeri- 
cal efficiency, a cubic approximation to a Gaussian, with 
standard deviations O.lg^/n and 0.3gy/n. The results 
with a Gaussian density of states have been compared to 
this cubic approximation, and no significant differences 
occur. 

In the condensed case, a third pair of lines are visi- 
ble, (labelled (c) in figure |2J). These correspond to the 
minimum energy for a neutral excitation, 2g^/nipQ , i.e. 
flipping the spin on a single site. They are analogous to a 
particle-hole excitation in BCS, which have twice the en- 
ergy of the gap for adding a single particle. This energy 
is the smallest "inter-band excitation" . If the inhomoge- 
neous broadening is small compared to the gap then these 
extra lines will be less visible, as seen in panel 2. For large 
inhomogeneous broadening , the peak at the edge of the 
gap (c) and the upper polariton (b) will merge, as shown 
by (e) in panel 4. This will result in a band of incoherent 
luminescence separated from the coherent emission at the 
chemical potential by the gap. Such structure, although 
associated with the internal structure of a polariton, can 



be seen even at densities where the transition tempera- 
ture is adequately described by a model of structureless 
polaritons. 

It is also interesting to consider the uncondensed cases. 
In panel 1, with the smaller broadening, three lines are 
visible; the upper and lower polariton, and between them 
luminescence from excitons weakly coupled to light (la- 
belled (a)). These weakly coupled states arise from exci- 
tons further away from resonance with the photon band. 
Although there may be a large density of such states, 
they make a much smaller contribution to luminescence 
than the polariton states, because of their small photon 
component. With larger broadening, these weakly cou- 
pled excitons form a continuum stretching between the 
two modes, as shown in panel 3. 



G. Momentum distribution of photons 

From the Green's function for photon fluctuations, one 
can calculate the momentum distribution of photons in 
the cavity. For a two-dimensional cavity coupled via 
the mirrors to three-dimensional photons outside, the 
momentum distribution may be observed experimentally 
from the angular distribution^. 

When uncondensed, N(p) is given by 



N(p) = lim{^( T + ^(T)) 

5 — >0+ 



/dx 
—n B (z)e Sz g 1 i(iz,p). 
2m 



(29) 



However when condensed, since the system is two- 
dimensional, it is necessary to treat fluctuations more 
carefully. Writing ip(r) — \/ po + n(r)e^( r \ the action 
involves only derivatives of the phase, showing that large 
phase fluctuations are possible. For quadratic fluctua- 
tions, one can use the matrix in eq. (|25|l describing trans- 
verse and longitudinal modes, and relate these to the 
phase and amplitude excitations. 

At low enough temperatures, it is possible to calcu- 
late N(p) by considering only the phase moda^i. Thus, 
neglecting amplitude fluctuations gives; 



N(p) = ± J d 2 r J d 2 r>e*i r - r ')po(e<^~^) 

,i J d 2 R J d 2 te ik.t e -D(R+t/2,R-t/2) ; (30) 



Po- 



where po is taken as the mean-field photon density. The 
phase correlator, found by inverting the amplitude/phase 
action, is: 



2itl3p a h l 



In 
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(31) 



The thermal length is £y = P c i where c is the veloc- 
ity of the sound mode from eq. I|23|) . This comes from 
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FIG. 2: Incoherent luminescence vs momentum (x axis) and 
energy (y axis) calculated above (panels 1 and 3) and below 
(panels 2 and 4) phase transition. Panels 1 and 2 are for inho- 
mogeneous broadening of O.lg^/n, and 3 and 4 are for 0.3g\/n. 
All other parameters are the same as in the unbroadened case 
shown in panel 1 of figure Q In order to map the infinite 
range of intensities to a finite scale, the colour is allocated as 
the hyperbolic tangent of intensity. Energies are in units of 
g^/n, and wavevectors in units of y/n. 



the energy scale at which fluctuations become cut by the 
thermal distribution. 

In this approximation, eq. (|30(1 may be evaluated 
exactly^ giving 



N(p) 



9 (p£r)" 
2vrp n — 



-2 1 " 



2-n 



P 



x 1 r/ Jo(x)dx 

„ r(i-y2) 

T(r,/2) ' 



(32) 



where rj = m/2Trf3poh controls the power law decay of 
correlations. The second line follows from an identity 
(see ref. exercise XVII. 32). This is valid only for 
small r], far from the transition. The Kosterlitz-Thouless 
transition^ occurs when rj becomes large, an approxi- 
mate estimate of the transition is at 77 = 2. 

The momentum distribution can therefore be calcu- 
lated both at low temperatures, where such a scheme 
holds, and at high temperatures, when uncondensed. 
This is shown in figure[3] When condensed, the power law 
divergence leads to a peak normal to the plane. This peak 
reflects coherence between distant parts of the system, so 
in a finite system this peak is cut at small momenta™. 
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FIG. 3: Momentum distribution of photons, from which fol- 
lows the angular distribution. These are plotted for low tem- 
perature condensed systems, where the approximation of in- 
cluding only phase fluctuations is valid, and for the simpler, 
uncondensed case. The inset illustrates the choices of density 
and temperature. (Parameters A* = I, m* = 0.01, wavevec- 
tor plotted in units of gy/n, temperature in units of g\/n, 
density in units of n, and N(k) in arbitrary units.) 



IV. FLUCTUATION CORRECTION TO THE 
MEAN-FIELD THEORY 

In this section, we calculate the fluctuation correc- 
tions to the mean-field density, and thus to the mean- 
field phase boundary. Our method is similar to that 
of Nozieres and Schmitt-Rink42, who studied fluctuation 
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corrections to the BCS mean-field theory for a model of 
interacting, propagating fermions. However, because our 
model differs from that of Nozieres and Schmitt-Rink, 
our approach to including second order fluctuations will 
also differ. We begin by presenting a brief summary of 
the method used by Nozieres and Schmitt-Rink, and fur- 
ther developed by Randeria 3 £. We then discuss how our 
approach differs from their work. 



in our model are of particular importance. Firstly, our 
boson field is dynamic, and there exists a chemical po- 
tential for bosons. In this respect, our model is closer to 
the boson-fermion models2£*Si studied in the context of 
Feshbach resonances, for which Ohashi & Griffin 2 — have 
studied fluctuation corrections. Secondly we consider a 
two-dimensional system; this requires calculation of fluc- 
tuations in the presence of a condensate, as discussed in 
the next section. 



1. Fluctuation corrections in three dimensions 

To consider fluctuation corrections to a mean field- 
theory, one first needs to find the partition function in 
terms of a coherent state path integral for a bosonic field, 



Z = J Vipexp(-S[i(j]) 



For the Dicke model, S[ip] is given in eq. JSJ. In the 
work of Nozieres and Schmitt-Rink, S[ip] resulted from 
decoupling a four-fermion interaction, and then integrat- 
ing over the fermions. This effective action may be un- 
derstood as a Ginzburg-Landau theory, with coefficients 
that are functions of temperature and chemical potential 
as well as the parameters in the Hamiltonian. To find 
the mean-field phase boundary, one needs both to find 
the values of temperature and chemical potential where 
the transition occurs, and to calculate the density evalu- 
ated at these parameters. 

For the mean-field theory, the action is evaluated for a 
static uniform field, an d minimised w.r.t. this field. 
At the critical temperature, a second order phase transi- 
tion occurs, and the minimum action moves to a non-zero 
■00- The density is found by differentiating the free en- 
ergy w.r.t. chemical potential. For the mean field theory, 
the free energy is approximated by the action evaluated 
at ipQ. 

To go beyond the mean field theory, one can expand 
the effective action about the saddle point: 



-S[V>o 



VSrJj exp — - — — - 



1 (PS_ 

2 dip 2 



(33) 



■00 



This gives an improved estimate of the free energy, from 
which follows an improved estimate of the density, and 
thus of the phase boundary. In three dimensions (but 
not in two, as discussed below), one only needs an esti- 
mate of the density at the mean-field critical tempera- 
ture. It is therefore sufficient to consider an expansion 
about the normal state saddle point, ip = 0. Such fluc- 
tuations may be understood as the contribution to the 
density from non-condensed pairs of particles, whereas 
the mean-field estimate of density included only unpaired 
fermions. These corrections will increase the density at 
a fixed critical temperature, or equivalently decrease the 
critical temperature for a given density. 

Because of features of our model, our approach differs 
from that of Nozieres and Schmitt-Rink; two differences 



2. Fluctuations in two dimensions 

To find the fluctuation-correction to the mean-field 
phase boundary in two dimensions, it is necessary to con- 
sider fluctuations in the presence of a condensate. Con- 
sidering fluctuations in the normal state would lead to 
the conclusion that the normal state can support any 
density: As one approaches the mean-field critical tem- 
perature, the fluctuation density will become infra-red 
divergent, allowing any density. This correctly indicates 
that no long-range order exists in two dimensions; how- 
ever a Kosterlitz-Thouless^Siii transition does occur. 

Instead one must start by considering fluctuations in 
the presence of a condensate. This gives a density, de- 
fined by the total derivative of free energy w.r.t. chemical 
potential, of the form: 



P : 



dF 



dF dF dipo 



(34) 



When considering fluctuation corrections in the pres- 
ence of a condensate, in any dimension, one must take 
care to consider the depletion of the order parameter due 
to the interaction between condensed and uncondensed 
particles. This is discussed in detail in section llV Al Such 
a depletion means that, for a fixed temperature, the crit- 
ical value of the chemical potential changes; at the mean- 
field critical chemical potential the formula for total den- 
sity may become negative. It is therefore necessary to 
make a separate estimate of the order parameter in the 
presence of fluctuations, and define the phase boundary 
where the order parameter goes to zero. 

In three dimensions such a calculation can be achieved 
by identifying parts of the density as the population of 
the ground state and of fluctuations (as discussed below 
in section IIV A 2 b|l . In two dimensions, no true con- 
densate exists, but a quasi-condensate with a cutoff fco 
can be considered. As discussed by Popo\*^, the quasi- 
condensate and fluctuation densities both contain terms 
which diverge logarithmically as fco — > 0, but these diver- 
gences cancel in the total density. 

Instead, in two dimensions, one must consider an al- 
ternate definition for the location of the phase bound- 
ary. Since the transition is a Kosterlitz-Thouless transi- 
tion, one should map the problem to the two-dimensional 
Coulomb Gas 43 . This requires the vortex core energy and 
strength of vortex-vortex interactions, which both scale 
as h 2 p s /2m, where p s is a superfluid response density. 



10 



The phase transition thus occurs when p s = #2mT/h 2 . 
The numerical pre-factor depends on the vortex core 
structure. However, approximating the critical condition 
by p s = leads to only a small shift to Tc- 



the fluctuation contribution to the second term in eq. I|34|) 
will take the form 



P = 



III 



\#o 



dp ' 



(38) 



A. Total derivatives and negative densities 

This section discusses the effects and interpretation of 
the second term in eq. H34JI . In ref. |23, Ohashi & Griffin 
define the density as the partial derivative of free energy 
w.r.t. chemical potential, neglecting the second term in 
eq. Ij34(l , which they describe as a higher-order correction 
under the Gaussian fluctuation approximation. As shown 
below in sec. II V ATI the contribution of the second term 
in eq. (|34|l to the density should not necessarily be ne- 
glected in the Gaussian fluctuation approximation. Sec- 
tion llV A2l shows explicitly that the two terms in eq. 134|) 
are of the same order. The existence of the second term 
is crucial in finding a finite density in two dimensions, 
however may be less important in three dimensions. 



Gaussian fluctuation approximation 



In section HvX2l we will show that the second term in 
eq. (|34|l is of the same order as the first. Before this, we 
explain why the second term of eq. I|34J) should not be 
automatically neglected. Even though it is of the form 
d 3 S/dijjQ, such terms are not necessarily small, and can 
contribute to the density at quadratic order. 

To see that a Gaussian theory may be correct even if 
such third order terms are not small, consider an expan- 
sion of the effective action, 



s = siM + S# 2 + S^ 3 + 

d% d% 



(35) 



A Gaussian approximation is justified if, using this ac- 
tion, the expectation of the cubic term is less than the 
quadratic term. This condition can be written as 



d 2 S 
d^ 



3/2 d 3 S 

> 

, / 1 1 



(36) 



This need not require that the coefficient of the cubic 
term is smaller, i.e. that 

d 2 S d 3 S 

In fact, if both terms are of the same order, but large; 
d 2 S/dtp 2 ~ d 3 S/dtjjQ > 1, then the condition is 
fulfilled. 

Writing the free energy including fluctuations schemat- 
ically as 



F = S[ip ] +ln 



XVtyexp ( ~^?<5V 2 



(37) 



where (. . .) signifies averaging over the fluctuation ac- 
tion. Thus, in calculating the condensate density, there 
is a term which depends on d 3 S/d^pQ but only involves 
second order expectations of the fields. Since d 3 S/dijj 3 ) 
is not necessarily small, and it contributes to the density 
at quadratic order, there is no a priori argument for ne- 
glecting this term. In the following we show explicitly 
that this term should be included. 



2. Total derivatives for a dilute Bose gas 

The following discussion will show explicitly that the 
two terms in equation (|34() are of the same order for a 
weakly interacting dilute Bose gas (W.I.D.B.G.), 

H-pN = ^(e fc - M)4 a fc + jE a i+ q a l'-q a k a k'- 

k k,k' ,q 

Further, the terms in eq. Ij34(l will be interpreted by con- 
sidering the Hugcnholtz-Pincs relation at one loop order, 
as discussed in ref. ^3 

a. Saddle point and fluctuations. To find the free 
energy, consider the static uniform saddle point, (aja ) = 
\A\ 2 = p/g, and quadratic fluctuations, which are gov- 
erned by the Hamiltonian: 

H eS = ^(e fc - p + 2gA 2 )a\a k 



9 A " (J t 



-k^ a k a -k 



(40) 



Thus, by a Bogoliubov transform, the free energy and 
density become 

F = - f iA 2 + - A 4 

+ E [\ ln (! - e ^ Ek ) + - ^ - M)) ,(41) 



A 2 -J2[n B (E, 



E k 2E k 



(42) 



where E k = ^e k {^k + 2p). 

The total density is thus less than the saddle point A 2 , 
and could be negative. From the form of the fluctuation 
Hamiltonian, eq. (|4(J|) . using gA 2 = p, it can be seen that 
the fluctuation contribution is: 

Pf = ~ E { ( a k a k) + \ (( a W-k) + (a k a^ k )) 



(43) 
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Using the results: 



a l a k 



It? \ e k + P . e k + p- E k 
n B {E k )—- 1 — , (44) 



2E k 



it is clear this matches ea. (|42ll . 

In contrast, taking partial derivatives, and neglecting 

the second term in eq. (fMj) gives Pf — J2k (^ a k a kj- I n 
two dimensions, for /i ^ 0, this expression will be infra- 
red divergent, while eq. H42() is not. 

b. Hugenholtz- Pines relation To identify the mean- 
ing of the terms in eq. 143|) . one can consider the 
Hugenholtz-Pincs relation for the normal and anomalous 
self energies A(u>,k), B(co,k) respectively: 



A(0,0)- 5(0,0) =fx, 



(46) 



The approximations in the previous section are equiv- 
alent to evaluating A and B at one loop order. As ex- 
plained by Popo\^, this becomes 

(i = 2g(p + pi) - g(p + Pi) 

- 2g 2 p Y J {G{k)G{-k)-g 1 {k)Q 1 {-k)). (47) 



one must identify the depleted condensate density, 
eq. I0SJ with 

p dip dFfi uct 

Pa = j KT, — ■ (51) 

g dp, oipo 

The derivatives w.r.t. the order parameter therefore 
describe a depletion of the order parameter due to fluc- 
tuations. Physically, interactions between the conden- 
sate and the finite population of non-condensed particles 
(at finite temperature) push up the chemical potential. 
With such a theory, there now exists a region of param- 
eter space where which is not condensed, but p > 0. In 
such a region it is essential to include modifications of 
the particle spectrum due to interactions to describe the 
normal state. 

c. Comparison of methods The phase boundary for 
the W.I.D.B.G. model with static interactions is pecu- 
liarly insensitive to the calculation scheme. This can be 
seen by consider the Hugenholtz-Pines relation at the 
transition. In general, the anomalous self energy van- 
ishes at the transition, so p — A(0, 0). Since po = 0, the 
total density is pi, which may be found from the fluctu- 
ation Green's function, p = ^ w k G(uj, k), 



G(w, k) = [iw + e fe - A(0, 0) + A(u, k)\ 



(52) 



Here po is the new condensate density, pi the density of 
non-condensed particles, and p\ is an anomalous particle 
density, pi — ^2 k Gi{k), with Qi the anomalous Green's 
function. The last term is a second order correction due 
to the three boson vertices of the form gA(aJa'a + a'aa). 
This last term in eq. I|47|) can be evaluated to be 2gp~i, 
leading to the result 



p . _ . 

Pa = (2pi +Pi), 
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(48) 



showing that po + p\ is less than the saddle point density. 

Compare this expression for the total density to that 
from saddle point and fluctuations, 



P = Po + Pi 



dp 



dp dipo 



(49) 



where p s , p . — p/g is the saddle point expectation of the 
density, and -Ffluct is the free energy from the fluctuation 
contributions. Since pi, the density of non-condensed 
particles can be identified as 



Pi = 



E 



dF au 



dp 



(50) 



where p = A(0, 0) has been used. If the self energy is 
static, A(cu, k) — A(0, 0), then at the transition the quasi 
particles are exactly free. Any approximation scheme 
which gives B(Q, 0) = when po = will then reproduce 
this result. For this reason, partial derivatives will give 
correct calculations of the phase boundary for a dilute 
Bose gas, but this does not remain true for dynamic self 
energies. 

For a boson-fermion model, such as polaritons, because 
of the dynamic self energy, total and partial derivatives 
will give different answers. For the three-dimensional 
case studied by Ohashi and Griffin, this will lead to crit- 
ical temperatures differing by a numerical factor, but in 
two dimensions using partial derivatives gives divergent 
answers. Were one to use partial derivatives, the density 
calculated from the condensed and non-condensed phases 
would agree at the critical chemical potential. However, 
for the total derivative, the density calculated at the new 
critical potential need not agree with that calculated from 
the non-condensed phase. A difference between these re- 
sults reflects the fact that both are approximations of 
the phase boundary, and is indicative of the Ginzburg 
criterion. 
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B. Total density for condensed polaritons 

From the effective action, eq. I|12f) . the free energy per unit area, including quadratic fluctuations may be written 



as: 



F = S[j> ] 1 f°° d 2 k 
A A (3 J (2tt)2 



(\ilu + hu k + K^uj)] 2 - \K 2 (u)\ 

i r d 2 k 



1 



d 2 k 



J (27T) 2 



In 



sinh(^!/2) sinh(/36/2) 
smh(/3E) smh{f3hu k /2) 

I 



(2ir) 2 



In 1 - 



a4E 2 



2 (to fc £: 2 - ftJ> e 2 ) 



(53) 



(54) 



Here E, a and £1,2 are defined as in section ITTTl The 
first term in braces is S[i/jo], the second and third to- 
gether are the fluctuation corrections. As discussed in 
section III Al the interaction has been cut off at a scale 
K m , so for k > K m , the action is that of a free gas of 



The total density is then given by: 



A 



1 - — tanh(/3E) 



K -d 2 k 
(2n) 2 



J 



photons, i.e. the second term in braces. In the third 
term, there are contributions both due to the Matsubara 
sum of eq. I|17fl . and due to the 5^ terms in eq. (|12|) . as 
discussed further in the appendix IA 21 



i "| r°° d 2 k 

f[^]-f[2E] + -+g(k)j + J k —n B (Juj k ), (55) 



where 



M = 

9(k) = 
C = 



1 



1 



dC 



2(3 (1 -C)dfj,' 

f3sech 2 (PE)9 2 n g 2 \A\ 2 
2 (HCj k E 2 - hu Q i 2 ) A ' 



In going from eq. (|54|) to eq. I|55ll the two integrals 
have been re-arranged, the second now describing only 
the free, high energy photons. Again, the last term, g(k), 
arises due to the 5 U terms. 

The derivatives of polariton energies that arise in cal- 
culating the density may be given in terms of the expres- 
sions A(k), B(k) as defined in eq. l(T%|) : 



26 



d/i. 



4£i 



dA{k) 
dfi 

2A{k) 



± 



1 



,/A(k) 2 - B{k) 
dA(k) dB(k) 



dA(k) 

dfi 
dB(k) 

dfi 



8E 
1G 



dE 
dfi 

h 2 k 2 



2m 



dfi 

2huj k - 4e 

dE 
2E—uj k - 
dfi 



dfi 
2%Cjq, 



ehCbn 



;2 



(56) 
(57) 
(.58) 



To find dE j dfi in the presence of a condensate, one can 



differentiate the gap equation, eq. ©, giving 



- 1 = 



dE 



g n 



dfi 2E 2 



(PEsech 2 ((3E) - tanh^-B)) . (59) 



C. Two dimensions, Superfluid response 

Having found an expression for the total density in- 
cluding fluctuations, it is necessary to consider how fluc- 
tuations change the critical chemical potential. As dis- 
cussed in the introduction to this section, in two di- 
mensions this requires consideration of the Kosterlitz- 
Thoulcss phase transition. The phase boundary is ap- 
proximated from the condensed state by the chemical 
potential at which the superfluid response vanishes. We 
therefore must calculate the normal and superfluid re- 
sponse in the presence of a condensate. 



1. Calculating normal response density 

Following the standard procedur o 44 i 45 , we consider the 
current, 



Tik 



JfoO) =J2~ ^k-q/2^k+q/2- 



(60) 



For a perturbation 8H = Ml. J, the linear response may 
be written (Ji(<?, 0)) = Xiji^^jW)- By symmetry, the 
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most general response function is, 



Xii(q) = XL(q)-^ 



+ XT(q) 8 



(61) 



By gauge symmetry, eq. <|6(J|) is a conserved current. It 
follows, via a Ward identity, that qtXijiq) — <Z.j'Tr(Cf (fc)), 
so the total density is p = mxL(q ^ 0)/A In contrast, 
the transverse response depends on only the density of 
normal particles, p n ormai = m-XT(q — > 0)/A 

In order to calculate the total density correctly, it is 
necessary to introduce vertex corrections, i.e. 



This form of the vertex correction means only the lon- 
gitudinal response is affected. Therefore, the standard 
procedure is to calculate the total density directly from 
the free energy, as in section llV Bl and the normal density 
by linear response. 

To one loop order, and neglecting vertex corrections, 
the response function is given by, 



Xtj(q->0) 



^^^(k^k)^), (64) 



X iM - o) = Tr (r ( (k, k)0(k) 7j (k, k)g(k)) 



(62) 



where 7»(p, q) = <J?,{pi + qi)/2m, and Ti(p,q) is chosen 
to satisfy the Ward identity. However, the diagrams re- 
quired to satisfy the Ward identity (see ref. Eih take the 
form, 



and so, taking the continuum limit, the normal density 
is given by 



Pn 



d k ft k 1 x — > , _ , _ , 

--^Tr(a(k)<r 3 a(k)<r 3 ). (65) 



(2tt) 2 2m (3 



Tj(p,q) = 7i(P>q) + (.Pi - Qi)f(p,q)- 



(63) 



J 



For the polariton system, the trace can be evaluated to give 



Pn 
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Co(fc) • (67) 



Again, in evaluating the Matsubara sum, one must consider the S u terms. The term Co(fc) is the difference between 
the true term at lo — 0, and the analytic continuation appearing in the Matsubara sum in eq. (|67|l . and is given by: 



C Q (fc) =2ax 
with a as defined in eq. (|15() . 
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(68) 



2. Total photon density 

For the polariton system there is an added complica- 
tion. Equation l|67|) gives the density of normal photons, 
but equation 155|) is the total excitation density (includ- 
ing excitons). It is therefore necessary to calculate the 
total photon density. 

This can be done by considering separate chemical po- 
tentials for photons and excitons, which are set equal 
at the end of the calculation. This means making the 
change, 



pN — ► p c 
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(69) 



This density is given by equation (|55|l with two 
changes. Firstly, the mean-field exciton density, 



1 - — tanh(/3£) 



(70) 



should be removed. Secondly, in f[x],g(k) derivatives 
should be taken w.r.t /i pn ot.- This means replacing equa- 
tions JsTJ) and (JSBJ> by: 



dA{k) 
dpphot. 
dB(k) 
dpphot. 



dE 

= 8E- 2hco k - 4?, 
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2 1,2 



= 16- 
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2m 



dE , 
2E—u k - E 2 
dp 



e 2 



(71) 
(72) 



in the action. The photon density is then total derivative 
w.r.t. the photon chemical potential, p p hot.- 



This makes use of the fact that d£ , /d/i p hot. = dE/dp, as 
can be seen from the gap equation. 
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V. PHASE BOUNDARY INCLUDING 
FLUCTUATIONS 

Combining the results of seetion llVl the phase bound- 
ary is found by plotting the total density (eq. lj55|> 1 at 
the value of chemical potential where the normal photon 
density (eq. matches the total photon density (dis- 

cussed in section llV C 2|) . The phase boundaries found in 
this way are plotted in figure 0] The form of the phase 
boundary can be explained by considering how, at finite 
temperatures, the occupation of excited states of the sys- 
tem depletes the condensate. Which excited states are 
relevant changes with density. 



Resonant case 




BECof 
photons^ ' 



BCS-like regime 



BEC of polaritons 

o.oi o'.i 
Density 



10 



When condensed, the lowest energy mode is the phase 
mode, described by eq. At low density, this has a 

shallow slope, and consequently a large density of states. 
Such excitations are described in a model of point bosons. 
The phase boundary can therefore be estimated from the 
degeneracy temperature of a gas of polaritons, of mass 
2m, where to is the bare photon mass: 



2irh up 



(73) 



As the density increases, the phase mode becomes 
steeper, and so has a smaller density of states. The rele- 
vant excitations are then single particle excitations across 
the gap. Such excitations are accounted for in the mean- 
field theory. Combining equations © and JJJ gives the 
result: 



T c = g^fn 



V 1 - 



2tanh~ 1 (l - 2p/n) - ln(p/n) ' 



(74) 



As seen in figure 0] the mean-field boundary is effec- 
tively constant on the scale of the boundary for BEC of 
point bosons, and so the crossover to mean-field always 
occurs near T ss g^/n, the Rabi splitting. The density at 
which this crossover occurs depends on the photon mass. 
Comparing equations (|75|l and ifT^) , this crossover occurs 

at a density p cross0 vor ~ «™*- _ 

In terms of the measurable Rabi splitting, g^/n and 
polariton mass to, this gives the density: 



mg^/n 



^crossover 



(75) 



Yamamoto et aZi 2 i 3 i 4 , 
TOeiectron- These values 



For the structures studied by 
g\fn w 7meV and to w 10~ 5 
give a crossover density of ^crossover = 2.6 x 10 8 cm~ 2 . 
This is both much less than the estimates of exper- 
imentally achieved density, n ~ 10 n cm~ 2 per pulse, 
and also much less than the Mott density in this struc- 
ture, nMott ~ 3.6 x 10 13 cm -2 . For the structures stud- 
ied by Dang et alm& g^/n ss 13meV, and to « 3 x 
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FIG. 4: Mean-field phase boundary, and phase boundaries in- 
cluding fluctuation correction for four values of photon mass, 
on a logarithmic scale (Insets are on a linear scale). The 
top panel is the resonant case, A = 0, adapted from ref. 
and shown for comparison. In the lower two panels the exci- 
ton is detuned below the photon band by A = 1.5gy/n and 
A = 2.5(?y / n respectively. Temperature plotted in units of 
g\fn and density in units of n. 
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10 _5 TO c i oc tron, so crossover densities are again of the same 
order, p C rossover = 5 x 10 8 cm~ 2 . Dang et al. also presents 
results for the detuned case 5 , discussed below, with a 
range dimensionless detunings 0.5 > A* > —0.7. 

Equation (|75l) describes the crossover in terms of prop- 
erties measurable for a given microcavity. However, to 
understand what fundamental lengthscales control this 
crossover density, it is necessary to write the coupling 
strength and polariton mass in terms of the dimensions 
of the cavity and properties of the excitons. Using the 
expressions in section [H] for photon mass and coupling 
strength g, this gives the crossover density as: 



crossover 



= 4tt 2 



1 d ah ^fn 



(76) 



The crossover density is therefore controlled by two 
parameters: The width of the cavity, w, and the ratio of 
electron-hole separation to average two-level system sep- 
aration, dab^/n = dab /r scp&ration- If the average two-level 
system separation (r sepa ration) is less than the electron- 
hole separation (d a fc), then our model of localised two- 
level systems will break down. Therefore, within our 
model the largest possible crossover density scale is 1/w 2 . 
This lengthscale occurs because the cavity size controls 
the wavelength of the lowest radiation mode. Crossover 
to a BCS-like mean-field regime occurs when the density 
approaches a scale set by the wavelength of light, rather 
than one set by the exciton Bohr radius. Therefore, in 
general this crossover density is much less than the Mott 
density. 

At yet higher densities, the single particle excitations 
are saturated, and so the condensate becomes photon 
dominated. In this regime, the transition temperature is 
that for a gas of massive photons. If the photon mass 
is large (to* > 1), a mean-field regime never exists, in- 
stead the phase boundary changes directly from polariton 
condensation to a photon condensation. However, since 
for experimental parameters the dimensionless mass is 
only of the order of 10 -3 , a mean-field regime will exist. 
These various crossovers are illustrated schematically in 
figure [SJ 



B. Effects of detuning 

If the excitons are detuned below the photon (positive 
detuning), it becomes possible for the system to reach 
half filling while remaining uncondensed. For positive de- 
tunings greater than 2gy/n the mean-field phase bound- 
ary becomes re-entrant, as shown in the bottom panel 
of figure 0] For smaller but still positive detunings, the 
mean-field boundary has a maximum critical density at 
a finite temperature, but no maximum of critical tem- 
perature. The opposite case, of excitons detuned above 
photons, shows no interesting features; the system will 
always condensed before half filling. 

This multi-valued phase boundary is discussed in 
ref. fill and can be explained in terms of phase locking 




FIG. 5: A schematic picture of how the relevant excitations 
change between the polariton BEC, the BCS-like mean-field 
and the photon BEC regimes. In the low density limit, a 
shallow sound mode exists. At higher densities, this becomes 
steeper, and the relevant excitations are gapped single par- 
ticle excitations. At yet higher densities, these modes are 
saturated, and the high k photon modes become relevant. 



of precessing spins*, either about spin down (low den- 
sity) or spin up (high density) states. Above inversion, 
increasing the density reduces the extent to which a spin 
may precess. For very large detunings, at low temper- 
atures, the phase diagram therefore becomes symmetric 
about half filling. 

When A > 2g^fn, the re-entrance leads to a point at 
zero temperature where two second order phase bound- 
aries meet. Including fluctuations, as shown in the bot- 
tom panel of figure ^ these phase boundaries no longer 
meet, but instead there is a region where two different 
condensed solutions coexist. In this region, there are two 
minima of the free energy, so we expect there will be 
a first-order phase boundary between them. Although 
this boundary could be calculated by comparing the free 
energies including fluctuations, its form may be altered 
significantly by higher order corrections. 

In the mean-field theory, at zero temperature, the 
chemical potential jumps discontinuously at the point 
where the two phase boundaries meet. This can be un- 
derstood by the chemical potential locking to the lower 
polariton for the lower density transitions, and to the 
upper polariton at higher densities. This can be seen 
in the lower panel of figure which plots the value of 
the chemical potential at the phase boundary. Including 
fluctuations, the jump in chemical potential has a simi- 
lar form, and is somewhat larger. In the region of coex- 
istence discussed above, the two minima of free energy 
have different chemical potentials, so at the first-order 
transition, the chemical potential will jump. 

At smaller detunings, as shown in the central panel of 
figure^ although the mean-field phase boundary is sin- 
gle valued, adding fluctuations can reproduce the same 
coexistence regions. For this to occur, the photon mass 
must be large — i.e. there must be a significant density 
of states for fluctuations. As shown in the upper panel 
of figure this coexistence is also characterised by two 
minima of the free energy, with different chemical po- 
tentials, and so is also expected to become a first-order 
transition in the same manner. 

To explain how fluctuations lead to the introduction 
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FIG. 6: Chemical potential vs density and temperature at 
the phase boundary, for the mean-field phase boundary, and 
m* = 0.50 fluctuation corrections. Plotted for detunings of 
A = l.bg^/n and A = 2.5g^/n, with all other parameters as 
in figure 0] Temperature and chemical potential plotted in 
units of g\/n, and density in units of n. 



of multiple phase boundaries at a single temperature, 
it is necessary to consider the upper branch of excita- 
tions, £,2{k). With positive detuning, the energy of this 
mode (w.r.t. chemical potential) can continue to fall as 
the chemical potential increases in the condensed state. 
This has two effects, it makes the sound velocity larger 
(as can be seen from eq. I)23f) '). and increases the popula- 
tion of this "pair-breaking" upper mode. The combina- 
tion of these effects is responsible for the creation of the 
coexistence region by fluctuations. 



C. Effects of inhomogeneous broadening 

It is interesting to consider how a small inhomogeneous 
broadening will modify the phase boundary. Exact cal- 
culations with a continuum of exciton energies are tech- 
nically challenging and not particularly illuminating, so 
the following presents a discussion of the main effects 
expected. The following discussion is for a Gaussian dis- 
tribution of energies, centred at the bottom of the photon 
band, with a standard deviation much less than g\fn. 



The most significant change to the boundary is due 
to the existence of a low energy tail of excitons. This 
means that, even at low densities, the chemical potential 
lies within the exciton band, and a BCS-likc form of T c 
will be recovered. Consider the density of states, 



v a {e) = 



cxp (-e 2 /2cr 2 



27TCT 



(77) 



For large negative chemical potentials, at low temper- 
atures and densities, the saddle point equation |(BJ| be- 
comes 
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(78) 



and the mean-field density Q , using the asymptotic form 
of the error function, is 



Pm.f. = v s \ — 



(79) 



The cutoff, A is approximately 2a 2 /p, but appears 
only as a pre-exponential factor, and so the density de- 
pendence it gives to T c will be neglected. Thus, the 
mean-field transition temperature at low densities then 
becomes 



T c = A cxp ( - 

9 P 



(80) 



For low densities, this result is very different to the 
mean- field theory without broadening, eq. I|74(l . Whereas 
before the mean-field boundary was approximately con- 
stant, it now drops rapidly to zero. If one now considers 
how fluctuations will modify this boundary, it is more 
helpful to consider the density as a function of temper- 
ature. At low temperatures, fluctuations increase the 
density by a small amount, Ap oc T. Without broad- 
ening, the mean-field critical density is approximately 
p w ne _9v ^/ T and goes to faster than the fluctua- 
tion corrections, Ap. Therefore, as shown in figure 01 
at low temperatures the fluctuation contribution con- 
trols the phase boundary. With broadening the mean- 
field critical density at low temperature is approximately 
p as 2cr 2 J g 1 ln(A/T), which goes to zero more slowly than 
Ap. Therefore, including fluctuations in this regime does 
not change the form of the phase boundary drastically. 
Hence at very low densities, the boundary is again well 
described by a mean-field theory. 



D. Relation to alternate models and experimental 
systems 

Our model describes two-level systems with a finite to- 
tal density of states — the density of states, per unit area, 
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integrated over all energies, is finite. As will be explained 
below, this finite total density of states is responsible for 
the re-entrant behaviour seen in the mean-field theory 
for detunings A > Ig^/n. In the preceding sections, the 
finite total density of states has also been implicated in 
explaining the existence of a photon dominated region at 
high temperature, and the multi-valued phase boundary 
in the presence of fluctuations. This section discusses 
how these effects may change in alternate models which 
do not have saturable two-level systems. Although the 
multi-valued phase boundary is expected only to occur 
for a finite band of two-level systems, the existence of a 
photon dominated regime at high densities is more gen- 
eral. 

Before discussing the more involved question of how 
changing the density of states affects fluctuations, we 
first summarise its effect on the mean-field theorjii 4 * 1 ^ 7 .. 
Consider the highest possible density achievable in the 
normal state. Since there exists a bosonic mode, the 
chemical potential cannot exceed the energy of this mode 
if the system is to remain normal. Therefore at zero tem- 
perature, only exciton modes below the boson mode are 
relevant. Regardless of whether the total density of states 
is finite, the density of states below the bosonic mode 
will be finite. However, at non-zero temperatures, exci- 
ton modes above the chemical potential can be occupied 
by the tail of the Fermi distribution. If there is a finite 
total density of states, there will be a maximum density 
of excitons that can be occupied thermally. This is what 
is referred to as "saturation" below. Note that for a ex- 
citon band centred at the bottom of the photon band, 
this maximum density of excitons is half the density of 
two-level systems — the system saturates at half filling. 
If the total density of states is not finite, it is possible to 
support any total density in the normal state by making 
the temperature high enough. 

This saturation is responsible for the multi-valued 
phase boundary in the mean-field theory. If the density 
is close to total inversion of the two-level systems, all 
two-level systems must be in the up state. This makes it 
hard for them to produce a macroscopic polarisation — 
viewed as spins, this is to say that if S z w +1/2, then 
S x will be small. Hence, the system becomes uncon- 
densed near total inversion. If the excitation density is 
greater than the density of two-level systems, the mean- 
field theory requires a coherent photon density. Hence, 
the system condenses again, giving a re-entrant bound- 
ary. Without saturable excitons, neither the uncondcns- 
ing due to reduction of mean-field polarisation, or the re- 
condensation due to exciton saturation will occur. With 
fluctuations, the multi-valued phase boundary is analo- 
gous to the mean-field re-entrance, and appears to require 
saturable excitations in the same manner. Therefore, we 
do not expect such multi-valued phase boundaries in a 
general model with a continuum of exciton states. 

Let us now consider the case where there is a finite den- 
sity of exciton states, separated by the exciton binding 
energy from a continuum of electron-hole states. In such 



a case, if the continuum is well separated from bound 
states, it may only affect the phase boundary at densi- 
ties larger than those where the features discussed above 
occur. Well separated here means that the exciton bind- 
ing energy is much larger than the energy scales in our 
model, in particular, much larger than g\fn. If the con- 
tinuum only has effects at very high densities, the exotic 
multi-valued phase boundary described in previous sec- 
tions will be realisable. For the systems studied by Ya- 
mamoto et al&, g\Jn ~ 7meV, and the exciton binding 
energy Ryd* lOmeV. For Dang et a/. 1 , gy/n w 13meV, 
and Ryd* w 25meV. In neither case can the effects of 
the continuum be avoided. Reducing the Rabi splitting, 
g^/n, might allow the multi- valued phase boundary to 
be observed. However, reducing the Rabi splitting will 
decrease the transition temperature in the region of inter- 
est. In addition, to have well strong coupling, the Rabi 
splitting must remain larger than the polariton linewidth 
due to photon lifetimes. 

The existence of a photon dominated region at high 
densities is however generic, and does not rely on a model 
with a finite total density of states. In an electron-hole 
plasma mode l 37 ' 48 , such a regime is also predicted. At 
large densities, as the chemical potential (lying within 
the exciton band) approaches the bottom of the photon 
band, the photon density increases much faster than the 
exciton density. The result in ref. |37J only describes a 
photon dominated regime at zero temperature. For the 
two-level system model, at high densities, the system re- 
mains photon dominated, and with increasing tempera- 
ture changes from coherent to incoherent photons. Such 
a change from coherent to incoherent photons is also ex- 
pected to occur in the electron-hole model, due to the 
large occupation of bosonic modes near the chemical po- 
tential, but further work is required here. 

This discussion of how our model relates to experimen- 
tal systems has so far concentrated on what happens at 
high densities. At lower densities (including the densities 
of current experiments), no such significant differences 
are expected. This is because, at low densities, higher en- 
ergy exciton modes would not be occupied, even if they 
exist. In particular, the angular distribution of radia- 
tion (figure |2J and excitation spectrum (figures ^ an d 0) 
should remain unchanged for generic models. Such sig- 
natures should therefore be expected for equilibrium con- 
densation in the systems currently studied. 

The signatures of condensation presented in this pa- 
per are calculated for a system in thermal equilib- 
rium, while current experiments are pumped. For non- 
resonantly pumped experimentsi^i&iiii, one must con- 
sider how pumping will modify the excitation spectrum 
and the occupation of modes. This can be described 
by coupling the system to baths describing pumping of 
excitons and decay of photons. For systems near equi- 
librium, with small coupling to baths, one expects the 
excitation spectrum to remain close to the equilibrium 
case, but with non-thermal occupations. Even with non- 
thermal occupation, the large density of states for exci- 
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tations near the chemical potential can be expected to 
produce a peak in the angular distribution of radiation. 
For strong coupling to baths, the spectrum of excitations 
will also change. One particularly significant change is 
the possibility of giving a finite lifetime to the Goldstonc 
mode. In this strongly pumped region, the signatures 
predicted in this paper can be expected to change signif- 
icantly. 

It is also of interest to discuss how these signatures 
are related to the behaviour seen in resonantly pumped 
cavities&^iifi. In these experiments pumping at a critical 
angle excites polaritons at momentum k p , causing emis- 
sion from signal, k = 0, and idler, k = 2k p , modes. Such 
a system may be described as an optical parametric os- 
cillator. Above a threshold, the luminescence from the 
signal increases superlinearly, and the signal linewidth 
narrows dramatically. In such a system, the occupation 
of the signal mode obeys a self-consistency condition, and 
the relative phase between the pump and signal mode is 
free. However, the nature of this self-consistency differs 
from that for an equilibrium condensate, and thus the 
signatures of equilibrium condensation are no longer im- 
mediately applicable. 

The form of the condensed luminescence spectrum, and 
the angular distribution of polaritons depend on the ex- 
istence of the low energy Goldstone mode. The energy of 
this mode vanishes as k — > as a consequence of the gap 
equation, eq. ©, which means that global phase rota- 
tions cost no energy. In a laser, the coherent field is also 
set by a self-consistency condition, balancing pumping 
and decay. Like condensation, the laser transition can 
also be described as symmetry breaking^. However, be- 
cause the self-consistency relates imaginary parts of the 
self energy, the dynamics of modes near the lasing mode 
is diffusive. Therefore a free global phase and a self- 
consistency condition are not sufficient for the signatures 
described in this paper. 

For the optical parametric oscillator experiments, the 
self-consistency equation is complicated by the existence 
of a coherent idler fiel d 50 i 51 . Since such experiments are 
strongly pumped it is expected that, despite the free 
phase and self-consistency, the luminescence spectrum 
and angular distribution of radiation as described in this 
paper will not be applicable. The laser and the equi- 
librium condensate are extreme cases, and the distinc- 
tion in practice is less clea»i2i^. For example, adding 
decoherence£2*£2i to an equilibrium condensate causes a 
crossover to a regime better described as a laser. 



VI. CONCLUSIONS 

We have studied the effect of fluctuations about the 
mean-field theory for a model of localised excitons cou- 
pled to a continuum of photon modes. When condensed, 
the presence of a gap in the fermion density of states, and 
the existence of the phase mode, cause dramatic changes 
to the spectrum of collective modes. Such changes lead 



to signatures of condensation in the luminescence and 
absorption spectrum (figure , and in the momentum 
distribution of radiation escaping the cavity (figure |3J) . 

Including the contribution to density due to fluctua- 
tions, we have studied the crossover from a BEC of po- 
laritons at low densities, through a BCS-like mean-field 
regime at intermediate densities, and finally to a BEC 
of massive photons at high densities, as shown in fig- 
ure The BCS-like regime occurs at densities achieved 
in current experiments. Our study of the crossover can 
be compared to other studies of the crossover; in Fesh- 
bach resonance systems2L£4, or for fermions interacting 
via a static four-particle interactionSiiii. In distinction 
to those systems, due to the nature of the fermion den- 
sity of states in our system, there is no clear difference 
between the roles of the number and gap equations in 
the BEC and BCS-like regimes. Rather, the crossover 
is in the nature of the fluctuations that depopulate the 
condensate. 

At low densities, fluctuation corrections significantly 
alter the form of the phase boundary from its mean-field 
form, leading to a dependence T cx p. As the density in- 
creases, the transition temperature approaches the Rabi 
splitting, and those single particle excitations which are 
included in the mean-field calculation dominate, leading 
to a recovery of the mean-field limit. This occurs at a 
density scale set by the wavelength of light, not the ex- 
citon separation, and so at densities much less than the 
Mott density. At yet higher densities, the system be- 
comes photon dominated and the boundary is that for a 
BEC of massive photons. If the exciton is detuned below 
the photon, the mean-field boundary can become multi- 
valued. With fluctuations, this multi-valued structure 
occurs for smaller detunings than are required for it to 
occur in the mean-field theory. At least for a saturable 
two-level systems, this multi- valued boundary likely in- 
dicates a first-order transition between two condensed 
states. 

Because our system is two-dimensional, it required 
fluctuations to be considered in the presence of the con- 
densate. In the presence of a condensate it is important 
to consider changes to the density both due to "conden- 
sate depletion" as well as the occupation of fluctuation 
modes. Such condensate depletion is included by taking 
full derivatives of the action w.r.t. the chemical potential. 
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APPENDIX A: LEHMANN REPRESENTATION 
AND BROKEN SYMMETRY 

1. Analytic properties of thermal and retarded 
greens functions 

The inverse thermal Green's function contains terms 
proportional to 5 U . Working from the Lehmann repre- 



sentation, it may be shown that such terms can arise 
in the thermodynamic Green's function, but not in the 
dynamic Green's functions. 



To see this, consider the standard Lehmann representation (see e.g. ref. |55|, section 17) for the retarded greens 
function for Bose fields: 

GrH = lim r P ^ X)d *, v Pl(x) = (l-e-^)^2\(n\i;\m)\ 2 e^ F - E ^S(x-E mn ). (Al) 

<5^0+ J_ OQ X - [U + id) ' 

However, for the thermal Green's function, an extra term may appear, 

g(u>) = [ dTe^T±(eP( F -^e HT ipe- HT tf)=Y f \{n\iP\m)\ 2 ^ (A2) 

J° V ' n,m J Q 

= f X ^^ + P6^\M\^\ 2 ^ F ~ En) ^n). (A3) 

J-oo X + IUJ f-^ 



The last term in l|A3(l can be identified as the contribu- 
tion to the Green's function due to transitions between 
degenerate states; or due to a macroscopic occupation of 
the photon in the ground state. Such a term does not 
occur for the retarded Green's function, and so in cal- 
culating the spectral Lehmann density pl{x), one may 
neglect its effects. 



2. Matsubara summation with thermal greens 
functions 



The 6 U terms will contribute to Matsubara sums in- 
volving the thermal greens functions, in properties such 
as the density. In performing the Matsubara summation 
one must use 
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where although A(z) involves /, the analytic continua- 
tion of /, /(0) does not involve analytic continuation. If 
the analytic continuation of / is regular at z = 0, then 



w " poles of / 

+ /(o)-/(o), 



/(=)fc„ th (| 



(A5) 



i.e. one must add a term to correct for the difference 
between / and its analytic continuation at z = 0. 
3. Emission and absorption coefficients 

Considering the Green's function for photon fluctua- 
tions; at zero temperature the function pl(x) would give 
the density of states, weighted by the photon component 
of a state. At finite temperatures, one may extract the 
probability to emit a photon, 

Pemit(x) = ^|(m|V|n)| 2 e^- E «)^ + S ron ) 

= n B (x)pi J (x), (A6) 

or to absorb a photon 

= (l + n B (x)) Ph (x), (A7) 



The energy x is measured w.r.t. the chemical potential, 
so at zero temperature there is only emission of photons 
at energies below the chemical potential, or absorption 
of photons above the chemical potential. The Lehmann 
density itself, ph(x) is the difference of these, and can 
be interpreted as an absorption coefficient, which when 
negative represents gain. 
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